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Abstract 

Nonlinear dynamics of a bouncing ball moving in gravitational field 
and colliding with a moving limiter is considered. Displacement of the 
limiter is a quadratic function of time. Several dynamical modes, such 
as fixed points, 2 - cycles and chaotic bands are studied analytically and 
numerically. It is shown that chaotic bands appear due to homoclinic 
structures created from unstable 2 - cycles in a corner-type bifurcation. 

1 Introduction 

Vibro-impacting systems belong to a very interesting and important class of 
nonsmooth and nonlinear dynamical systems [TJ f2J [31 S] with important tech- 
nological applications [5J [HI [Zl IE] • Dynamics of such systems can be extremely 
complicated due to velocity discontinuity arising upon impacts. A very charac- 
teristic feature of such systems is the presence of nonstandard bifurcations such 
as border-collisions and grazing impacts which often lead to complex chaotic 
motions. 

The Poincare map, describing evolution from an impact to the next impact, 
is a natural tool to study vibro-impacting systems. The main difficulty with 
investigating impacting systems is in finding instant of the next impact what 
typically involves solving a nonlinear equation. However, the problem can be 
simplified in the case of a bouncing ball dynamics assuming a special motion 
of the limiter. Bouncing ball models have been extensively studied, see [9] and 
references therein. As a motivation that inspired this work, we mention study 
of physics and transport of granular matter [6] . A similar model has been also 
used to describe the motion of railway bogies [7] . 
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Recently, we have considered several models of motion of a material point 
in a gravitational field colliding with a limiter moving with piecewise constant 
velocity |1Q[ ITT1 [T2l IT51 114) . Moreover, we have proposed more realistic yet still 
simple models approximating sinusoidal motion of the table as exactly as possi- 
ble but still preserving possibility of analytical compuations 15 . In the present 
work we study the model in which displacement of the table is a quadratic and 
periodic function of time. 

The paper is organized as follows. In Section 2 a one dimensional dynamics 
of a ball moving in a gravitational field and colliding with a table is reviewed and 
the corresponding Poincare map is constructed. A bifurcation diagram is com- 
puted for displacement of the table assumed as quadratic and periodic function 
of time. In the next Section dynamical modes shown in the bifurcation diagram 
such as fixed points, 2 - cycles and chaotic bands are studied analytically and 
numerically. It is shown that chaotic bands appear due to homoclinic structures 
created from unstable 2 - cycles in a corner-type bifurcation. We summarize 
our results in Section 4. 



2 Bouncing ball: a simple motion of the table 

We consider a motion of a small ball moving vertically in a gravitational field 
and colliding with a moving table, representing unilateral constraints. The ball 
is treated as a material point while the limiter's mass is assumed so large that 
its motion is not affected at impacts. A motion of the ball between impacts is 
described by the Newton's law of motion: 

mx = —mg, (1) 

where x — dx/dt and motion of the limiter is: 

y = y(t), (2) 

with a known function y. We shall also assume that y is a continuous function 
of time. Impacts are modeled as follows: 

= y(n), (3) 

i(r+)~y(n) = -R(x(r-)-y(n)), (4) 

where duration of an impact is neglected with respect to time of motion between 
impacts. In Eqs. (j3)), ((4]) stands for time of the i-th impact while x~ , i^are 
left-sided and right-sided limits of Xi (t) for t ~ > Tj , respectively, and R is the 
coefficient of restitution, < R < 1 [5]. 

Solving Eq. |T|) and applying impact conditions dU we derive the 
Poincare map [16) : 

jY(T i+1 ) = 7 y(r i )-A? +1 +A i+1 ^ > (5a) 
Vi+i = -RVi+2RA i+1 + 7 (l + R)Y(T i+1 ), (5b) 
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where Aj+i = 7V|_i — T; t . The limiter's motion has been typically assumed 
in form Y S (T) = sin(T), cf. [11] and references therein. This choice leads to 
serious difficulties in solving the hrst of Eqs.([S]) for Tj+i, thus making analytical 
investigations of dynamics hardly possible. Accordingly, we have decided to 
simplify the limiter's periodic motion to make (|5ap solvable. 

In our previous papers we have assumed displacement of the table as piece- 
wise linear periodic function of time |11[ 1121 114) . In our recent work preliminary 
results for function Y(T) assumed as quadratic Y q and two cubic functions of 
time, Y Cl and Y C2 have been obtained [TH] . In this work we study dynamics for 
quadratic function of time Y q (T): 



y q {T) 



-WTlT 



o < f < | 



16 ( T - | 



T - 1 



< T < 1 



(6) 



with T = T — [T\ , where [x\ is the floor function - the largest integer less than 
or equal to x, in more detail. 

Since the period of motion of the limiter is equal to one, the map ([5]) is 
invariant under the translation Tj — > T< + 1. Accordingly, all impact times Tj 
can be reduced to the unit interval [0, 1]. The model consists thus of equations 
|5]), ([6]) with control parameters R, 7. 

In Fig. 1 below we show the bifurcation diagram with times of impacts 
computed for growing 7 and R = 0.85 (see also |15] where bifurcation diagram 
with velocities just after impacts against 7 was shown) . 
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Figure 1: Bifurcation diagram. R = 0.85, 7 G [0.0101, 0.0635]. 



It follows that dynamical system ([5]), ([6]) has several attractors: six fixed 
points, one 2 - cycle and, possibly, chaotic attractor. 
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3 Analytical results 

We shall investigate now fixed points, 2 - cycle and chaotic bands shown in Fig. 
1, combining analytical and numerical approach. 

3.1 Fixed points 

We shall first study periodic solutions with one impact per k periods. Such 
states have to fulfill the following conditions: 

V n+1 = V n = K (fc/1) , T n+1 = T n + k = Ti fe/1) +k (k = 1, 2, . . .) , (7) 

where: 

Ti fc/1) e (0, 1) , Vj k/1) > iY q (T) . (8) 
Substituting these conditions into ([5]), © we obtain two sets of fixed points: 

< T (s) = i - < I (q) 

U - 1 * ~ 4 32 7 (l + i?) - 2 V 3 ) 

V* = k 

where the impact occurs in time interval £ (0, |) an d 

1 <- rp{u) _ 3 , k{l-R) s 1 / ln \ 

2 - 1 * — 4 + 32 7 (l + fi) - 1 

V* = k 

with impacts taking place in time interval £ (h, l). 

Solutions (0) fulfill physical requirements and are stable in the following 
interval of 7: 

k 8(l+% —1— 8(1+V' 
(s) 

where lower bound is a consequence of T* > while the upper bound follows 
from the condition that eigenvalues A of the stability matrix obey |A| < 1. On 
the other hand, solutions (TTOj) are always unstable and are physical for: 

7>fc^y, (12) 
(u) 

what is equivalent to the condition T„ < 1. 

3.2 2 - cycle 

It follows from the bifurcation diagram, Fig. 1, that there exists a stable 2 
- cycle with time of first impact T*i g (0, ^) and time of the second impact 
T*2 G l). Such periodic solution must fulfill the following equations which 
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(13) 



are easily obtained from Eqs. ([5]), ©: 

f I67 (T.2 - |) (T* 2 - 1) = 

-167T»i (T»i - I) - (T* 2 - T*i) 2 + (T 1 ^ - T*i) Ki 
K2 = -RV*\ + 2R (T„ a - T»i) + 7 (1 + i?) (32T, 2 - 24) 
-167(T,3-1)((T,3-1)-^) = 

I67 (T* 2 - \) (T„2 - 1) - (T»3 - T* 2 ) 2 + (T* 3 - T* 2 ) K2 
K 3 = - W, 2 + 2i? (T.3 - T» 2 ) + 7 (1 + .R) (8 - 32 (T* 3 - 1)) 
T* 3 =T ml + l 
. K 3 = Ki 

Eliminating variables we arrive at equation for time of the first impact only: 
C4X 4 + C 3 a; 3 + C 2 x 2 + C x x + C = 0, (14) 

where a; = Xi* and 

^ /7a L11D2 o , -A (247 + 1287 2 + l)i? 3 + (3847 2 +87-l)i? 2 + (3847 2 -247 + l)fl+128 7 2 -87-l 

Co = ((87 + 1) R* - 87 + 1) i '- > (fl+ i)V.-i) 2 

(24 7 +1287 2 + l)ii: 4 -27?, 3 + (-5127 2 +2+327)fl 2 + (-2-5127 2 +647)_R-1287 2 +87+l 
Cl - " b4 ^ (fl-l^fl+l) 2 

C 2 = -2048 7 2 - fl3+( - 1 +S^-; +487)fl+1 

C 3 ^4096 7 2 (-l + 16 7 ) fl ( 2 +^ 1 
C 4 = -4096 7 2 (-1 + I67) f±i 

(15) 

It follows from the bifurcation diagram that the 2 - cycle is born when 
Ti = x = 0. This in turn occurs when Co = 0. Equation Co = has three 
roots: 

'1 8 1-R 2 ' 

_ 1 -3-R 2 +2fl+l-Vfl 4 -12H a -2_R 2 +4TR+9 f-i R \ 
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32(i?+l) R+1 

1 -3J? 2 +2fl+l + VJ? 4 -12_R a -2fl 2 +4_R+9 
32(i?+l) fi+1 

Testing Eqs. (| 16[) against numerical computations we find out that the stable 2 
- cycle is born at 7 = 73. For example, for R = 0.85 we have 73 = 0.023 3674, 
cf. Fig 1. 



3.3 Chaotic bands and homoclinic structure 

Magnification of the bifurcation diagram near the origin of chaotic bands is 
shown in Fig. 2 below. 

There are six chaotic bands (and six basins of attraction) above the critical 
point 7cr . The first band which appears for appropriate initial conditions is 
shown in Fig. 3. Each band consists of two subbands only since due to cyclic 
periodic conditions points T = 0, T = 1 are identified. 

We note that there is a switch of stability - fixed points become unstable 
precisely at 7 = 7 cr when chaotic bands appear. Then at 7 = "f cr a homoclinic 
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Figure 2: Bifurcation diagram. Six chaotic bands, R = 0.85, 7 € 
[0.0627, 0.0630]. 



trajectory with T = 1~ is probably created - this is suggested by presence 
of clusters of points near T = 0, T = 1. Computer simulations show that 
near 7 cr there are six unstable 2 - cycles. For 7 < -f cr impacts occur in the 
following time intervals T tl e (0, ±), T», 2 e (l, T* 3 G (2 + m, 2 + m + §), 
m = 0, 1, 2, . . ., and T»3 = T*i mod 1 where true impact times without cyclic 
conditions are shown. Then at 7 = 7 cr a corner event occurs, i.e. T*2 = 1. 
Indeed, at T»2 = 1 acceleration of the table is discontinuous. Then, for 7 > ^ cr 
another six unstable 2 - cycles are created with impact times T*i € (0, 
T* 2 € (5, 1), r*3 € (2 + m, 2 + m+ |), to = 0.1.2. . . ., and T* 3 = T rl mod 1. 

It is now possible to write down equations, suggested by numerical compu- 
tations, for the unstable 2 - cycles for 7 > 7 cr : 

I67 (y m ~ 5) (Vm - 1) = -167a; (x - |) - (y - a;) 2 + (y - a;) u 
v = + 2R (y - x) + 7 (1 + J?) (32t/ ro - 24) 
( -167Z„ (zn - |) = I67 (y m - |) (y m - 1) - (z - y m ) 2 + (z - y m )v ^ 
w = -Rv +2R(z- y m ) + 7 (1 + iZ) (8 - 32z„) ' ' 

^71 — X 

w = u 

where 

x = T»i, y = T* 2 , z = T« 3 , 
Um = T *2 - m, z n = T» 3 - n, 
u = V*i, w = V*2, w = Via, 

with integer to, ri where numerical computations suggest that n 



(18) 

= m + 2. 
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Figure 3: Bifurcation diagram. First chaotic band, R = 0.85, 7 € 
[0.0627, 0.0630]. 



We are going to solve equations ([P?]) , (fl"8]) at 7 = 7 cr and this means that we 
have to put y m = 1 . It follows that there are six unstable 2 - cycles in question 
which are obtained for m — 0, 1, . . . , 5 and n = m + 2, solutions for larger to's 
being nonphysical. Solving these equations we get: 

7~ = fW' " = ^ + 2, (19) 

^(cr) _ _ 2_R 2 -6-4m+2 A /(3+5m+2m 2 )fl 4 + l+3m+2m~ _ , . 

m = ^ ~ 2(J? 2 -1) — TO — z!, l^ZUJ 

and we do not show more complicated expressions for , , ■ 

It follows from Eq. (f2T)]) that for R = 0.85 there are only six acceptable 
solutions with T*i > corresponding to six chaotic bands in Fig. 2. In the 
Table 1 impact times and the corresponding velocities just after the impact, 
computed from Eqs. fT7)l. ([15]). (12TJ1) . are listed for i? = 0.85 and to = 0, 1, . . . , 5: 



Table 1 



TO 





1 


2 


3 


4 


5 


j,(cr) 
*1. m 


0.434 7 


0.343 6 


0.260 2 


0.178 4 


0.0974 


0.0167 


y(cr) 


0.5148 


1.623 7 


2.716 9 


3.806 5 


4.894 6 


5.9819 


- 1 *2, m 


1~ 


1- 


1" 


1- 


1" 


1- 


y{cr) 
v *2, m 


1.454 6 


2.366 7 


3.279 4 


4.192 3 


5.1051 


6.0180 



Critical 2 - cycle to = 0, n = 2 (R = 0.85) is shown below. 
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Figure 4: Critical 2 - cycle. Position of the 2 - cycle is denoted by open circles. 

In Figs. 5, 6 the second chaotic band and first two chaotic bands are shown. 
The bifurcation diagrams for 7 > j cr = 0.06291 . . . were computed for initial 
conditions shown in the Table 1. 

Sharp edges of chaotic bands are given within good approximation by T^ 1 "^ 
and also T = 0, 1. It seems that the homoclinic structure exists for all values of 
7 > 7cr shown in the Figures. 
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Figure 5: Bifurcation diagram. Second chaotic band, R — 0.85, 7 € 
[0.0627, 0.0630]. 
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Figure 6: Bifurcation diagram. First two chaotic bands, R - 
[0.0627, 0.0630]. 

We have also computed Lyapunov exponents, cf. the Table 2. 
Table 2 



0.85, 7 e 



m 





1 


2 


3 


4 


5 




0.2 


1.2 


1.6 


1.85 


2.0 


2.1 



4 Discussion and closing remarks 

We have studied dynamics of a material point moving vertically in a gravita- 
tional field and colliding with a limiter. Displacement of the limiter has been 
assumed as quadratic function of time ^ . Due to the simplicity of the problem 
it was possible to investigate the dynamics analytically with some support from 
numerical computations. Firstly, fixed points were found and their stability was 
determined. Secondly, equations for a stable 2 - cycle were found and simplified, 
cf. Eqs. (fi"4|) . (fl~5|) . From these equations analytical condition for birth of the 
2 - cycle was found (cf. 7 = 73 in Eq. tfTB"]) ). Finally, a transition to chaotic dy- 
namics was described in analytical terms. It was shown that six stable chaotic 
bands appear from six unstable 2 - cycles. Equations for these 2 - cycles were 
found and solved to yield critical value of 7, Eq. (1191) . and impact times and 
the corresponding velocities at 7 = 7 cr , see Eq. (|20| where T^ r m was given. 
Approximation to the band edges was also found. 

We have demonstrated, combining analytical and numerical approach, that 
at the transition point 7 = j cr unstable 2 - cycles give rise to homoclinic struc- 
tures which lead to chaotic behaviour. This transition is a corner-type bifurca- 
tion similar to that found in a bouncing ball model with piecewise linear velocity 
[I"4] . In our future work we shall study models with displacement of the table 
described by a cubic functions of time [15] . 
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